Poynting's theorem for planes waves at an interface: a scattering matrix approach 
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We apply the Poynting theorem to the scattering of monochromatic electromagnetic planes waves 
with normal incidence to the interface of two different media. We write this energy conservation 
theorem to introduce a natural definition of the scattering matrix S. For the dielectric-dielectric 
interface the balance equation lead us to the energy flux conservation which express one of the 
properties of 5": it is a unitary matrix. For the dielectric-conductor interface the scattering matrix, 
that we denote by S, is no longer unitary due to the presence of losses at the conductor. However, 
the dissipative term appearing in the Poynting theorem can be interpreted as a single absorbing 
mode at the conductor such that a whole S, satisfying flux conservation and containing S and this 
absorbing mode, can be defined. This is a simplest version of a model introduced in the current 
literature to describe losses in more complex systems. 

PACS numbers: 41.20.Jb, 42.25.Bs, 42. 25. Fx, 42.25. Gy 



I. INTRODUCTION 

The scattering matrix S" is a useful tool to describe 
the multiple scattering that occur when planes waves en- 
ter to a system which in general is of complex nature, 
like atomic nucleus, chaotic and/or disordered systems, 
etcjii^ Although it is well known for particle waves in 
quantum mechanics^i^i^ii^ it can also be applied to any 
kind of plane waves. In the electromagnetic context the 
transfer matrix^'^ is known instead of S'Z but they are 
equivalent and related^. By definition, S relates the out- 
going plane waves amplitudes to the incoming ones to the 
system, from which the reflection R and transmission 
T coefficients are obtained; they are called reflectance 
and transmitance in electromagnetismj^ In the absence of 
dissipation, as happens in quantum electronic and some 
electromagnetic systems, S becomes a unitary matrix, in 
particular R + T — 1. However, if the system contains a 
dissipative medium, S is no longer unitary, in fact it is a 
sub- unitary matrix and R+T < 1, where the lack of unity 
is called absorbance in the electromagentic subject.— 

We are concerned with the electromagnetic case. 
There, a natural definition of S through one of its prop- 
erties arises from the Poynting theorem. This theorem is 
a energy balance equation which in the simplest form, i. 
e. for linear and non dispersive media, is given by^ 



V-Sp(r,i) 



du{r, t) 
dt 



-ReJ(r,i) • ReE(r,i), (1) 



where 



Sp(r,t) =ReE(r,i) x ReH(r,t) 



(2) 



is the Poynting vector which gives the energy flux per 
unit area per unit time, and u(r, t) is the electromagnetic 
energy density 

M(r,t) = i[ReE(r,t) • ReD(r, t) + ReB(r, • ReH(r,t)]. 

(3) 

Here, we have assumed that the electric E(r,t), elec- 
tric displacement D(r, t), magnetic induction B(r, t), and 



magnetic H(r,f) fields, as well as the current density 
J(r, i), are complex vectors whose real part only has 
physical meaning. We have also written explicitely the 
dependence on the position r and time t. The term in 
the right hand side of Eq. ([1]) is the negative of the work 
done by the fields per unit volume and represents the 
conversion of electromagnetic energy to thermal (or me- 
chanical) energy. The version of Poynting's theorem for 
dispersive media will not be touched here. 



Our purpose in this paper is to apply the Poynting 
theorem to the simplest scattering system, an interface 
between two different media, to illustrate the relation of 
Poynting's theorem and one property of S. First, we 
will consider the absence of dissipation for the dielectric- 
dielectric interface for which the defined S matrix is 
unitary. Second, the dissipative case is considered in 
the dielectric-conductor interface for which the scattering 
matrix, called S, is sub-unitary. The last system is the 
simplest example to explain quantitatively a model, that 
we call the "parasitic channels" model, introduced in con- 
temporary physics^i^ to describe complex systems with 
losses.— Also, S could represents the scattering of an 
"absorbing patch" used to describe surface absorptionjii 



The paper is organized as follows. In the next sec- 
tion we write the time averaged Poynting's theorem and 
calculate the corresponding Poynting's vector for planes 
waves and its flux through an open surface in a dielec- 
tric medium, as to be used in the sections that follow 
Sect, ini The dielectric-dielectric interface is considered 
in Sect. IIIII while Sect. |lV]is devoted to the dielectric- 
conductor interface. Finally, we conclude in Sect. |V] 
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II. POYNTING'S THEOREM FOR LINEAR 
AND NON-DISPERSIVE MEDIA 

In what follows we will consider monocromatic high 
frequency oscillating fields such that 



E(r,i) = E(r)e^'''^* 
H(r,i) = H(r)e-*'^*, 
J(r,t) = J(r)e-^-', 



(4) 
(5) 
(6) 



where have written the precise dependence on the spa- 
tial and temporal variables to avoid confusion due to the 
abuse of the notation. Here, the temporal average is of 
importance. Using the definition 



(7) 

for the average of a time dependent function f{t) over a 
period r, the average of Eq. ^ can be written as 

V • (Sp(r,t)) = -(ReJ(r,t) • ReE(r,t)), (8) 

where we have used that, with help of Eqs. ^ and ([5]), 

r-r g 



dt 



u{r,t)dt = u{r,T)-u{r,0) ^0. (9) 



Also, by substitution of Eqs. dH) and (O into Eq. (g]), it 
is easy to see that 

(Sp(r,i)) =ReSp(r), (10) 

with (we write again the precise dependence) 



Sp(r) = -E(r) X H*(r), 



(11) 



being the time averaged Poynting's vector. In equivalent 
way, using Eqs. Q and we get 

(ReJ(r,t) • ReE(r,i)) = ^Re [J(r) • E*(r)] . (12) 

Finally, the time averaged energy flux conservation law 
is given by 



V.Sp(r) = --J(r).E*(r), 



(13) 



where only the real part has physical meaning. 

If we integrate over a volume V enclosed by a close 
surface S, it says that the net flux $ of Sp(r) through S 
is the negative of the work done by the fields if there are 
dissipative components in V. I. e., 



where, from one side, 



W^^ [ J(r) •E*(r)dF 
2 Jv 



(14) 



(15) 



and, for the other side, using the divergence theorem $ 
can be written as a surface integral over E as 



$ = ^ Sp(r) • nda. 



(16) 




S(0) 



nda 



FIG. 1: A plane wave with linear polarization in the x di- 
rection traveling in the z direction cross the open surface of 
area A. The long arrows represent the corresponding time 
averaged Poynting's vector. 



A. Flux of Poynting's vector for plane waves 
through an open surface of area A 

For a monocromatic plane wave with linear polariza- 
tion in the x axis, propagating along the positive z- 
direction, the spatial component is 

E(z) = Eo e'^"- X, with k = — (17) 

c 

being the wave number, n the index of refraction of the 
medium and c the speed of light in vacuum; x is a unit 
vector pointing in the positive cc-axis, and the complex 
number Eq is the amplitude of E(z). The magnetic field 
H(2;), calculated from 'Ei{z) using one of the Maxwell 
equations (Faraday's law of induction), is 



H(z) = E^e 



(18) 



where y is a unit vector pointing in the positive y- 
direction and /io is the permeability of the free space 
(we have assumed a non- magnetic medium). 

Sustituting Eqs. ^ and HH]) into Eq. ^ we get 



S(?)(z) = ^|i?0pz, 



(19) 



where z is a unit vector pointing in the positive z-axis. 
This equation means that the time averaged energy flux 
is constant along the propagation. If we consider an open 
surface of area A as in Fig. n the flux of S^°^ (z) through 



it is (see Eq. ([16] 



$0 = 



— I^ol'rfa- — l^oPA (20) 
where we have taken the normal unit vector of A as z. 



III. POYNTING'S THEOREM AT A 
DIELECTRIC-DIELECTRIC INTERFACE 

Let us consider an interface between two dielectrics, 
with refractive indices n and n', in the xy plane as shown 
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FIG. 2: Interface between two dielectrics with indices of re- 
fraction n and n' . The amphtudes for the incoming {Ea and 
Eai) and outgoing {Ei, and Ey) planes waves are shown with 
the corresponding Poynting's vectors represented by long ar- 
rows. Short arrows denote the unit vectors fii and £12 of the 
covers of the cylindrical surface E. The cross section of the 
cylinder is A. 



in Fig. [21 For the shake of simphcity, we consider plane 
waves with linear polarization in the x axis with normal 
incidence on both sides of the interface. Therefore, the 
spatial part of the electric fields, E{z) on the left and 
E'(z) on the right, are of the form 



E(z) 
E'(z) = (^Ea 



rp ^tkz I rp —ikz 



Ak' z 



Eh' e 



X 

ik' z 



(21) 
(22) 



where Ea and Ea' {Eb and Ebi) are the complex ampli- 
tudes of the incoming (outgoing) planes waves. 

For this case, the right hand side of Eqs. and 
is zero {W = 0) because there is not dissipation. Hence, 
$ = 0. For convenience, we take a cylinder of cross 
section A, including the covers, to be the closed surface 
E (see Fig. [2]). Due to normal incidence only the flux 
through the covers of section A contribute to <&. Then, 
Eq. HH) gives 



0, 



(23) 



where (to = a, a', b, b') is the flux through one of the 
covers with normal unit vector hj (j — I, 2) due to the 
TO-th plane wave; it is given by 



= / S(^")(z)-n,da. 
Ja 



(24) 



With the convention that hj points outwards of E, $f, 
and $f,' are positive quantities while $a and $a' are im- 
plicitely negative. Using the result of Eq. (pO| for a plane 
wave, Eq. (j23p can be written as 

n{\Eb\^^\Ea\')+n'{\Eb'\^-\Ea'\^)^0, (25) 
which can be arranged in a matricial form that we will 



use later, namely 

( Et El. ) 
= i E* E* 



n 
n' 



Eb 
Eb' 



n 
n' 



Ea 
Ea' 



(26) 



A. The scattering matrix S for the interface 

By definition the Fresnel coefficients relate the outgo- 
ing to the incoming plane waves amplitudes as^i^ii^ 



Eb 
Eb' 



Ea 
Ea' 



where Sp the 2x2 matrix 



-rp t'p 



(27) 



(28) 



The Fresnel coefficients of reflection r^, r'p and trans- 
mission tp, t'p are^ 



rp 
tp 



n + rt 

2n 
n' + n ' 



t'^ 



-rp 

In' 

n' + n 



(29) 
(30) 



Although for this particular case they are real numbers, 
in general they are complex, in which case Sp is a com- 
plex matrix. In order to be more general, in what follows 
we will assume that is complex. 

By substitution of Eq. ([77]) , Eq. (pS)) can be written as 



K K' ) 4 



K K' 



n 
n' 



n 
n' 



Ea 
Ea' 

Ea 
Ea' 



from which we see that 



Si 





^ 


) Sp = 


' n 


M 




n' 


[ 


n' I 



(31) 



(32) 



(33) 



(34) 



S is known in the literature as the scattering matrix which 
has the following general structure^ 



or, equivalently, 

S^S = I, 

where / is the 2x2 identity matrix and 
5 = 




S = 



r t' 
t r' 



(35) 
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Here, r and t {r' and t') are the reflection and transmis- 
sion amplitudes for incidence on the left (right). From 
Eqs. (IMl), (ESI), dSQl) and ^ it is easy to see that 



t = 



n — n 
n + n 



t' = t. 



n 



n' 



(36) 
(37) 



Two remarks are worthy of mention.^ The first one is that 
t' — t because the system described by S is invariant 
under inversion of time, but in a more general context 
that is not necessary the case.— The second is that r' ^ r 
because the optical path is not the same for the reflected 
trajectory when incidence is from right or left. Of course, 
Eq. is satisfied in the case we are considering as can 
be easily checked using Eqs. ([M]) and (|37p. In particular, 

i? + T=l, (38) 

where R = |rp and T = |tp are the reflection and trans- 
mission coefficients. 

Eqs. P7|) and (IM)) implies that by definition S relates 
the outgoing to incoming plane waves amplitudes but 
normalized with the index of refraction. I. e.. 



£b' 



S 



Ea- 



rn 



where £i = ^JnlEi {i = a, b, a' , b'), and the electric field 
on both sides is given by 



E(z) = ( ^ e''^- + A e-^kz 
In \ n 



E'(z) 



£a,' ^ik'z 

n' 



n' 



(40) 
(41) 



We recall that Eq. ([39]) is a definition of S that has 
arised in a natural way from the restriction imposed by 
flux conservation, Eq. (|33p . and its structure, given by 
Eq. (|35p . reflects the symmetries present in the problem, 
one of them being the time reversal invariance. 



IV. THE DIELECTRIC-CONDUCTOR 
INTERFACE 

When one of the two media is a conductor with electric 
conductivity a, the one on the right in Fig. [3] let say, the 
treatment is equivalent as in Sect. Illll but with a complex 
index of refraction: n' n' + ir], where n' and rj are the 
optical constants!^ Also, the corresponding wave number 
becomes complex: we replace k' k' -\- in with 



k' = n'bj / c and k = rjuj /c. 



(42) 



Then, on the conductor side the electric field is an evanes- 
cent wave while on the dielectric side we have incoming 
and outgoing plane waves. They are 

ikz I t;^ „~ikz\ 



E(z) - {Ea^ 



Ebe- 



E'(z) = Eb,e-'''e'''''ic. 



(43) 
(44) 



n +ir\ 




At 
X 



A 

z 



FIG. 3: The medium on the right is a conductor whose con- 
ductivity is a. The refraction index, as well as the wave vec- 
tor, becomes complex such that only evanescent waves are 
present with an amplitude decaying exponentially as z 



mcreases. 



From Eq. or Eq. dM]), with Ea' = we see that 

Eb ^ rEa (45) 

Eb' = J^^tEa, (46) 
y n' + irj 

where r and t are obtained from Eqs. ([55]) and ([57]) but 
replacing n' by n' + 177, namely 



IT] 



t 



n + n' + ifj 
2y^n{n' + ir]) 
n + n! -\- irj 



(47) 
(48) 



In this case, the scattering matrix of the system is 1 x 1 
and we denote it by S. By definition (it is not necessary 
to normalize with respect to the index n) 



Eb — S Ea- 

Comparing with Eq. ([i5|) we see that 

S^r = ^e'\ 
where, from Eqs. ([T71) . 

[n — n'Y -1-77^ —2nr\ 



R = 



tan( 



77^ 



(49) 
(50) 

(51) 



As for S*, the Poynting theorem impose a restriction to 
S. To apply Eq. (|14p we consider that the closed surface 
E in Fig.[3]extends to infinity on the right of the interface, 
such that Eq. gives 



o\EvYe 



a 

4^' 



(52) 



where we used that 3{z) = o'E'(z). Then, Eq. ([H]) can 
be written as 



n 



{\Eb\' 



\E„, 



a 

2k. 



Eb.? 



(53) 
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where the area A has been cancelled. Using Eqs. ([46 
and ((^ . after an arrangement, Eq. ([55)1 gives 



(54) 



meaning that S* is a 1 x 1 subunitary matrix with Tp the 
lack of unitarity {strength of absorption or absorbance), 



lK\fr. 



\t? = 



Ann! 



77^ 



[n + rt')2 + r\ 



(55) 



The last equality is valid for metals in the upper infrared 
part of the spectrum and for metals at microwave and 
lower frequencies!^ Of course Tj, = 1 — |5p as can be 
easily verified using Eqs. l(5T|) . 

Eq. ([51)1 can be seen as resulting of the unitarity condi- 
tion for an S'-matrix that satisfy flux conservation (com- 
pare with Eq. ((38|) ). This S'-matrix should has the struc- 
ture (see Eq. 



S = 



/Tpe 



Tpe 



1(4,-6) 



(56) 



where the phase can be taken as the phase of 
^/nt/y^n' + ir] (see Eqs. (gSl) and (|55l) ): using Eq. (|48p 



tan( 



V 



n' + n 



(57) 



This form of S says that the losses because of the con- 
ductor can be interpreted as due to a single mode of ab- 
sorption whose "coupling" to the interface is Tp. This is 
the simplest version of what we call the "parasitic chan- 
nels" model^i^ appeared in the literature of contemporary 
physics in recent years to describe power losses in more 
complex systems (see Ref. 10 and references there in). 
That model consists in simulate losses with Np equivalent 
absorbing modes, each one having an imperfect coupling 



Tp to the system. The total absorption is quantified by 
7 — NpTp but Np and Tp can not be determined sepa- 
rately. Our result not only explain this abstract model 
but quantify in an exact way the coupling Tp of each ab- 
sorbing mode as well the scattering matrix of a single 
absorbing patch in the surface absoprtion model. ^"'^ Our 
treatment presented here can also be used to construct 
artificially a system with multiple absorbing modes^i^ 



V. CONCLUSIONES 

We reduced the time averaged Poynting theorem to a 
property of the scattering matrix. For that we applied 
this balance equation to a simplest scattering system, 
consisting of normally incident planes waves at an inter- 
face between two media. The simplest version of this the- 
orem was used such that dispersive media and the corre- 
sponding dissipation were ignored. Two kind of interfaces 
were considered. In the first one, a dielectric-dielectric in- 
terface, the energy fiux conservation leads to a natural 
definition of the scattering matrix S which is restricted 
to be a unitary matrix. We recalled that the structure of 
S reflects the symmetries present on the problem in other 
contexts. In the second one, the dielectric-conductor in- 
terface, the definiton of S'-matrix was used to describe 
the scattering taking into account the losses on the con- 
ductor side via the Poynting's theorem. This allowed us 
introduce the parasitic channels model used in contem- 
porary physics to describe the scattering with losses in 
more complex systems. We should were able to quantify 
the coupling of the single parasitic mode of absorption 
in our simple system. A system with multiple absorbing 
modes can be constructed and the results will be pub- 
lished elsewhere. Finally, the same treatement can be 
generalized for oblique incidence. 
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